Chronic Myeloid Leukemia: Molecular Monitoring of Residual Disease by Genomic DNA Compared to Conventional mRNA Analysis in Follow-Ups up to 8 Years

Chronic myeloid leukemia (CML) is a clonal myeloproliferative disorder resulting from the t(9;22)(q34;q11) balanced reciprocal translocation within a pluripotent stem cell (SC). The resulting Philadelphia (Ph) chromosome produces BCR-ABL1 fusion gene coding for a deregulated Abl tyrosine- kinase with constitutive and tumorigenic activity. The first line therapy of CML is imatinib mesylate, which targets Bcr-Abl protein, inhibiting proliferation pathways. Complete cytogenetic response can be achieved in 95% of patients treated in the early chronic phase (CP)1. Molecular monitoring of minimal residual disease is crucial to detect poor responses to imatinib and optimizing treatment with second generation tyrosine-kinase inhibitors or allogeneic stem cell transplantation. Residual leukemia is assessed by a quantitative manner evaluating levels of BCR-ABL1 transcripts by real-time reverse transcriptase PCR (qRT-PCR). Although qRT-PCR detects mRNA levels in a very sensitive manner, the negative result is difficult to interpret, because undetectable levels of chimeric transcript can reflect either an effective elimination of leukemia cells, or the presence of a quiescent leukemia SC transcriptionally silent

BCR-ABL1 genomic DNA PCR response kinetics during first-line imatinib treatment of chronic myeloid leukemia

Accurate quantification of minimal residual disease during treatment of chronic myeloid leukaemia guides clinical decisions. The conventional minimal residual disease method, RQ-PCR for BCR-ABL1 mRNA, reflects a composite of the number of circulating leukemic cells and the BCR-ABL1 transcripts per cell. BCR-ABL1 genomic DNA only reflects leukemic cell number. We used both methods in parallel to determine the relative contribution of the leukemic cell number to molecular response. BCR-ABL1 DNA PCR and RQ-PCR were monitored up to 24 months in 516 paired samples from 59 newly-diagnosed patients treated with first-line imatinib in the TIDEL-II study. In the first 3 months of treatment BCR-ABL1 mRNA values declined more rapidly than DNA. By 6 months the two measures aligned closely. The expression of BCR-ABL1 mRNA was normalized to cell number to generate an expression ratio. The expression of e13a2 BCR-ABL1 was lower than that of e14a2 transcripts at multiple time points during treatment. BCR-ABL1 DNA was quantifiable in 48% of samples with undetectable BCR-ABL1 mRNA, resulting in minimal residual disease being quantifiable for an additional 5-18 months (median 12 months). These parallel studies show for the first time that the rapid decline in BCR-ABL1 mRNA over the first 3 months of treatment is due to a reduction in both cell number and transcript level per cell, whereas beyond 3 months falling levels of BCR-ABL1 mRNA are predominantly due to depletion of leukaemic cells.Ilaria S. Pagani, Phuong Dang, Ivar O. Kommers, Jarrad M. Goyne, Mario Nicola, Verity A. Saunders, Jodi Braley, Deborah L. White, David T. Yeung, Susan Branford, Timothy P. Hughes, and David M. Ros

Germline mutations in mitochondrial complex I reveal genetic and targetable vulnerability in IDH1-mutant acute myeloid leukaemia

The interaction of germline variation and somatic cancer driver mutations is underinvestigated. Here we describe the genomic mitochondrial landscape in adult acute myeloid leukaemia (AML) and show that rare variants affecting the nuclear- and mitochondriallyencoded complex I genes show near-mutual exclusivity with somatic driver mutations affecting isocitrate dehydrogenase 1 (IDH1), but not IDH2 suggesting a unique epistatic relationship. Whereas AML cells with rare complex I variants or mutations in IDH1 or IDH2 all display attenuated mitochondrial respiration, heightened sensitivity to complex I inhibitors including the clinical-grade inhibitor, IACS-010759, is observed only for IDH1-mutant AML. Furthermore, IDH1 mutant blasts that are resistant to the IDH1-mutant inhibitor, ivosidenib, retain sensitivity to complex I inhibition. We propose that the IDH1 mutation limits the flexibility for citrate utilization in the presence of impaired complex I activity to a degree that is not apparent in IDH2 mutant cells, exposing a mutation-specific metabolic vulnerability. This reduced metabolic plasticity explains the epistatic relationship between the germline complex I variants and oncogenic IDH1 mutation underscoring the utility of genomic data in revealing metabolic vulnerabilities with implications for therapy.Mahmoud A. Bassal, Saumya E. Samaraweera, Kelly Lim, Brooks A. Bernard, Sheree Bailey, Satinder Kaur, Paul Leo, John Toubia, Chloe Thompson-Peach, Tran Nguyen, Kyaw Ze Ya Maung, Debora A. Casolari, Diana G. Iarossi, Ilaria S. Pagani, Jason Powell, Stuart Pitson, Siria Natera, Ute Roessner, Ian D. Lewis, Anna L. Brown, Daniel G. Tenen, Nirmal Robinson, David M. Ross, Ravindra Majeti, Thomas J. Gonda, Daniel Thomas, Richard J. D, Andre

Impact of additional genetic abnormalities at diagnosis of chronic myeloid leukemia for first-line imatinib-treated patients receiving proactive treatment intervention

Early view: March 23, 2023The BCR::ABL1 gene fusion initiates chronic myeloid leukemia (CML), however evidence has accumulated from studies of highly selected cohorts that variants in other cancer-related genes are associated with treatment failure. Nevertheless, the true incidence and impact of additional genetic abnormalities (AGAs) at diagnosis of chronic phase (CP)-CML is unknown. We sought to determine whether AGAs at diagnosis in a consecutive imatinib-treated cohort of 210 patients enrolled in the TIDEL-II trial influenced outcome despite a highly proactive treatment intervention strategy. Survival outcomes including overall survival, progression-free survival, failure-free survival and BCR::ABL1 kinase domain mutation acquisition were evaluated. Molecular outcomes were measured at a central laboratory and included major molecular response (MMR, BCR::ABL1 ≤0.1%IS), MR4 (BCR::ABL1 ≤0.01%IS) and MR4.5 (BCR::ABL1 ≤0.0032%IS). AGAs included variants in known cancer genes and novel rearrangements involving the formation of the Philadelphia chromosome. Clinical outcomes and molecular response were assessed based on the genetic profile and other baseline factors. AGAs were identified in 31% of patients. Potentially pathogenic variants in cancer-related genes were detected in 16% of patients at diagnosis (including gene fusions and deletions) and structural rearrangements involving the Philadelphia chromosome (Ph-associated rearrangements), detected in 18%. Multivariable analysis demonstrated that the combined genetic abnormalities plus the ELTS clinical risk score were independent predictors of lower molecular response rates and higher treatment failure. Despite a highly proactive treatment intervention strategy, first-line imatinib-treated patients with AGAs had poorer response rates. This data provides evidence for the incorporation of genomically-based risk assessment for CML.Naranie Shanmuganathan, Carol Wadham, NurHezrin Shahrin, Jinghua Feng, Daniel Thomson, Paul Wang, Verity Saunders, Chung Hoow Kok, Rob M. King, Rosalie R. Kenyon, Ming Lin, Ilaria S. Pagani, David M. Ross, Agnes S.M. Yong, Andrew P. Grigg, Anthony K. Mills, Anthony P. Schwarer, Jodi Braley, Haley Altamura, David T. Yeung, Hamish S. Scott, Andreas W. Schreiber, Timothy P. Hughes and Susan Branfor

Geographic population structure analysis of worldwide human populations infers their biogeographical origins

The search for a method that utilizes biological information to predict humans’ place of origin has occupied scientists for millennia. Over the past four decades, scientists have employed genetic data in an effort to achieve this goal but with limited success. While biogeographical algorithms using next-generation sequencing data have achieved an accuracy of 700 km in Europe, they were inaccurate elsewhere. Here we describe the Geographic Population Structure (GPS) algorithm and demonstrate its accuracy with three data sets using 40,000–130,000 SNPs. GPS placed 83% of worldwide individuals in their country of origin. Applied to over 200 Sardinians villagers, GPS placed a quarter of them in their villages and most of the rest within 50 km of their villages. GPS’s accuracy and power to infer the biogeography of worldwide individuals down to their country or, in some cases, village, of origin, underscores the promise of admixture-based methods for biogeography and has ramifications for genetic ancestry testing

An interlaboratory study of TEX86 and BIT analysis of sediments, extracts and standard mixtures.

Two commonly used proxies based on the distribution of glycerol dialkyl glycerol tetraethers (GDGTs) are the TEX86 (TetraEther indeX of 86 carbon atoms) paleothermometer for sea surface temperature reconstructions and the BIT (Branched Isoprenoid Tetraether) index for reconstructing soil organic matter input to the ocean. An initial round-robin study of two sediment extracts, in which 15 laboratories participated, showed relatively consistent TEX86 values (reproducibility ±3-4°C when translated to temperature) but a large spread in BIT measurements (reproducibility ±0.41 on a scale of 0-1). Here we report results of a second round-robin study with 35 laboratories in which three sediments, one sediment extract, and two mixtures of pure, isolated GDGTs were analyzed. The results for TEX86 and BIT index showed improvement compared to the previous round-robin study. The reproducibility, indicating interlaboratory variation, of TEX86 values ranged from 1.3 to 3.0°C when translated to temperature. These results are similar to those of other temperature proxies used in paleoceanography. Comparison of the results obtained from one of the three sediments showed that TEX86 and BIT indices are not significantly affected by interlaboratory differences in sediment extraction techniques. BIT values of the sediments and extracts were at the extremes of the index with values close to 0 or 1, and showed good reproducibility (ranging from 0.013 to 0.042). However, the measured BIT values for the two GDGT mixtures, with known molar ratios of crenarchaeol and branched GDGTs, had intermediate BIT values and showed poor reproducibility and a large overestimation of the "true" (i.e., molar-based) BIT index. The latter is likely due to, among other factors, the higher mass spectrometric response of branched GDGTs compared to crenarchaeol, which also varies among mass spectrometers. Correction for this different mass spectrometric response showed a considerable improvement in the reproducibility of BIT index measurements among laboratories, as well as a substantially improved estimation of molar-based BIT values. This suggests that standard mixtures should be used in order to obtain consistent, and molar-based, BIT values

Multi-wavelength observations of blazar AO 0235+164 in the 2008-2009 flaring state

The blazar AO 0235+164 (z = 0.94) has been one of the most active objects observed by Fermi Large Area Telescope (LAT) since its launch in Summer 2008. In addition to the continuous coverage by Fermi, contemporaneous observations were carried out from the radio to γ-ray bands between 2008 September and 2009 February. In this paper, we summarize the rich multi-wavelength data collected during the campaign (including F-GAMMA, GASP-WEBT, Kanata, OVRO, RXTE, SMARTS, Swift, and other instruments), examine the cross-correlation between the light curves measured in the different energy bands, and interpret the resulting spectral energy distributions in the context of well-known blazar emission models. We find that the γ-ray activity is well correlated with a series of near-IR/optical flares, accompanied by an increase in the optical polarization degree. On the other hand, the X-ray light curve shows a distinct 20 day high state of unusually soft spectrum, which does not match the extrapolation of the optical/UV synchrotron spectrum. We tentatively interpret this feature as the bulk Compton emission by cold electrons contained in the jet, which requires an accretion disk corona with an effective covering factor of 19% at a distance of 100 R g. We model the broadband spectra with a leptonic model with external radiation dominated by the infrared emission from the dusty torus. © 2012. The American Astronomical Society. All rights reserved

Insights into the high-energy γ-ray emission of Markarian 501 from extensive multifrequency observations in the Fermi era

We report on the γ-ray activity of the blazar Mrk 501 during the first 480 days of Fermi operation. We find that the average Large Area Telescope (LAT) γ-ray spectrum of Mrk 501 can be well described by a single power-law function with a photon index of 1.78 ± 0.03. While we observe relatively mild flux variations with the Fermi-LAT (within less than a factor of two), we detect remarkable spectral variability where the hardest observed spectral index within the LAT energy range is 1.52 ± 0.14, and the softest one is 2.51 ± 0.20. These unexpected spectral changes do not correlate with the measured flux variations above 0.3 GeV. In this paper, we also present the first results from the 4.5 month long multifrequency campaign (2009 March 15-August 1) on Mrk 501, which included the Very Long Baseline Array (VLBA), Swift, RXTE, MAGIC, and VERITAS, the F-GAMMA, GASP-WEBT, and other collaborations and instruments which provided excellent temporal and energy coverage of the source throughout the entire campaign. The extensive radio to TeV data set from this campaign provides us with the most detailed spectral energy distribution yet collected for this source during its relatively low activity. The average spectral energy distribution of Mrk 501 is well described by the standard one-zone synchrotron self-Compton (SSC) model. In the framework of this model, we find that the dominant emission region is characterized by a size ≲0.1 pc (comparable within a factor of few to the size of the partially resolved VLBA core at 15-43 GHz), and that the total jet power (≃1044 erg s-1) constitutes only a small fraction (∼10-3) of the Eddington luminosity. The energy distribution of the freshly accelerated radiating electrons required to fit the time-averaged data has a broken power-law form in the energy range 0.3 GeV-10 TeV, with spectral indices 2.2 and 2.7 below and above the break energy of 20 GeV. We argue that such a form is consistent with a scenario in which the bulk of the energy dissipation within the dominant emission zone of Mrk 501 is due to relativistic, proton-mediated shocks. We find that the ultrarelativistic electrons and mildly relativistic protons within the blazar zone, if comparable in number, are in approximate energy equipartition, with their energy dominating the jet magnetic field energy by about two orders of magnitude. © 2011. The American Astronomical Society

A method for next-generation sequencing of paired diagnostic and remission Samples to detect mitochondrial DNA mutations associated with leukemia

Somatic mitochondrial DNA (mtDNA) mutations have been identified in many human cancers, including leukemia. To identify somatic mutations, it is necessary to have a control tissue from the same individual for comparison. When patients with leukemia achieve remission, the remission peripheral blood may be a suitable and easily accessible control tissue, but this approach has not previously been applied to the study of mtDNA mutations. We have developed and validated a next-generation sequencing approach for the identification of leukemia-associated mtDNA mutations in 26 chronic myeloid leukemia patients at diagnosis using either nonhematopoietic or remission blood samples as the control. The entire mt genome was amplified by long-range PCR and sequenced using Illumina technology. Variant caller software was used to detect mtDNA somatic mutations, and an empirically determined threshold of 2% was applied to minimize false-positive results because of sequencing errors. Mutations were called against both nonhematopoietic and remission controls: the overall concordance between the two approaches was 81% (73/90 mutations). Some discordant results were because of the presence of somatic mutations in remission samples, because of either minimal residual disease or nonleukemic hematopoietic clones. This method could be applied to study somatic mtDNA mutations in leukemia patients who achieve minimal residual disease, and in patients with nonhematopoietic cancers who have a matched uninvolved tissue available.Ilaria S. Pagani, Chung H. Kok, Verity A. Saunders, Mark B. Van der Hoek, Susan L. Heatley, Anthony P. Schwarer, Christopher N. Hahn, Timothy P. Hughes, Deborah L. White and David M. Ros

Early BCR-ABL1 kinetics are predictive of subsequent achievement of treatment-free remission in chronic myeloid leukemia

With treatment-free remission (TFR) rapidly becoming the ultimate goal of therapy in chronic myeloid leukemia (CML), there is a need to develop strategies to maximise sustained TFR by improving our understanding of its key determinants. Chronic phase CML patients attempting TFR were evaluated to identify the impact of multiple variables on the probability of sustained TFR. Early molecular response dynamics were included as a predictive variable, assessed by calculating the patient-specific halving time of BCR-ABL1 after commencing tyrosine kinase inhibitor (TKI) therapy. Overall, 115 patients attempted TFR and had ≥12 months follow-up. The probability of sustained TFR, defined as remaining in major molecular response off TKI therapy for 12 months, was 55%. The time taken for the BCR-ABL1 value to halve was the strongest independent predictor of sustained TFR: 80% in patients with a halving time of 21.85 days (last quartile) (P<.001). The e14a2 BCR-ABL1 transcript type and duration of TKI exposure before attempting TFR were also independent predictors of sustained TFR. However, the BCR-ABL1 value measured at 3 months of TKI was not an independent predictor of sustained TFR. A more rapid initial BCR-ABL1 decline after commencing TKI also correlated with an increased likelihood of achieving TFR eligibility. The association between sustained TFR and the time taken for BCR-ABL1 to halve after commencing TKI was validated using an independent dataset. These data support the critical importance of the initial kinetics of BCR-ABL1 decline for long-term outcomes.Naranie Shanmuganathan, Ilaria S. Pagani, David M. Ross, Sahee Park, Agnes S.M. Yong, Jodi A. Braley, Haley K. Altamura, Devendra K. Hiwase, David T. Yeung, Dong-Wook Kim, Susan Branford, Timothy P. Hughe